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Abstract 

We employ the method of the theory of open quantum systems to analyze 
spin relaxation and decoherence in semiconductors in the presence of a mag- 
netic field. We derive a set of Bloch equations for electron spin with a fully 
microscopic determination of longitudinal and transverse relaxation times. 
Electron scattering from optical and acoustic phonons and random impurities 
is taken into account. We obtain explicit expressions for the spin relaxation 
times in terms of material constants and coupling strengths, exhibiting for- 
mal agreement with earlier treatments in the zero magnetic field limit with 

microscopic specification of their phenomenological parameters. 
PACS numbers: 03.65.Yz, 72.25.Rb, 72.15.Lh 
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I. INTRODUCTION 



Modern statistical physics encompasses a wide range of scientific problems from elemen- 
tary particles to the evolution of the Universe. Its methods facilitate the determination of 
the macroscopic characteristics of complex nonequilibrium systems without having to ex- 
actly analyze the highly complicated individual microscopic dynamics of each constituent 
part. In many cases the system of interest can be separated into two distinct subsystems. 
One of them, which we will term the " dynamical subsystem" , has only few degrees of free- 
dom, whereas the other, the " bath" , has , in an ideal case, the infinite number of them. The 
bath is a source of random perturbations (fluctuations) for the dynamical subsystem, also 
functioning as an absorber of dissipative energy. This separation is at the core of the theory 
of open quantum systems (TOQS) based on papers of Schwinger |]], developed in Refs. @-|4| 
and successfully applied to radiation damping || and some problems in solid state physics 

In the present paper we demonstrate the applicability of TOQS to the problem of spin 
dynamics in semiconductors, which has been extensively studied over the last decade both 
experimentally and theoretically || in conjunction with numerous proposals for a vari- 
ety of spin-based devices ||. The most promising materials for device purposes, the III-V 
and II- VI compounds, have been shown to have spin relaxation rates dominated by the 



D'yakonov-Perel' (DP) mechanism fH| at moderate temperatures and low hole concentra- 
tions. In contrast to earlier treatments, we develop a fully microscopic stochastic theory 
of spin dynamics in the presence of an external magnetic field, taking account of pertinent 
scattering mechanisms. Our analysis is based on an innovative implementation of the two 
step relaxation process corresponding to the relaxation time hierarchy involved in (a) elec- 
tron thermalization, and (b) spin relaxation. In the first stage of solution, we determine 
the relaxation rates and fluctuation characteristics of electron orbital motion due to cou- 
pling to optical and acoustic phonons and random impurities as a bath. Spin relaxation 
dynamics (the slowest process in the system) can be neglected in stage (a). The second 
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stage, (b), proceeds with the analysis of the spin relaxation process due to spin-orbit inter- 
action, wherein the orbital degrees of freedom are considered as an "effective heat bath", 
having the characteristics determined in stage (a). The orbital motion can be considered as 
an intermediary transmitting fluctuations from the phonon/impurity bath to electron spin, 
and, also, transferring the dissipated energy flow in the opposite direction. A set of Bloch 
equations with two distinct relaxation times (longitudinal relaxation time, T±, responsible 
for spin magnetic moment relaxation, and transverse relaxation time, T 2 , responsible for 
decoherence) is derived in this second stage. In both stages of our analysis we employ the 
TOQS method. 

The structure of the paper is as follows. In Section II we present a brief outline of TOQS 
and the principal relations used in subsequent Sections. The microscopic formulation of the 
spin relaxation problem is given in Section III. In Section IV we analyze the orbital electron 
dynamics in the presence of an external magnetic field with coupling to optic and acoustic 
phonons and random impurities. The spin-orbit interaction is examined in Section V, where 
we derive the Bloch equations on a microscopic basis. A summary of this work is given in 
Section VI. 

II. THEORY OF OPEN QUANTUM SYSTEMS 

The full Hamiltonian of the system, separated into a "dynamical subsystem" and a 
"bath", can be written as 



where H is the Hamiltonian of the dynamical subsystem, Hb is the Hamiltonian of the 
bath and H int describes their interaction. We can write the last term in a product form 



where F a (t) is, in the general case, a nonlinear function of the dynamical subsystem variables, 
Q a (t) is a function of the bath variables, and the summation index a can refer to either 



H = Hq + Hb + Hi 



(1) 




(2) 



a 
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coordinate projections or a mode index. One can see from Hi nt that both the dynamical 
subsystem and the bath influence each other. Moreover, F a plays a role of a generalized 
force conjugate to the generalized bath coordinate, Q a and vice versa (Figurel). 

In the case of Gaussian statistics of the unperturbed bath variables, or when the coupling 
of the dynamical subsystem to the bath is weak, the full Heisenberg operator of the bath is 
given by f| 



is a linear response function, or the retarded Green's function of the unperturbed bath 
variables. r\ (t — ti) is the Heaviside unit step function. Another important function is the 
correlation function of the unperturbed bath variables 



Here and below we use the notation [...,...]_ for a commutator, and [...,...], for an anticom- 
mutator. The linear response function of Eq.(|J), and the correlation function of Eq.(||]), are 
related by means of the fluctuation-dissipation theorem 




(3) 



where 




(4) 




(5) 




B 



(6) 



where 




(7) 



and 





To derive the fluctuation-dissipation theorem we employed the relation 
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(IQI(0),Q>)} + ) = WMQMUcoth^, (9) 

based on the Gibbs distribution and, therefore, the fluctuation-dissipation theorem, in the 
present form, is valid only for systems in equilibrium. 

The Heisenberg equation of motion for an arbitrary operator of the dynamical subsystem, 
A(t), is given by 

Mt) = ^l A (t),H }- - i[A(t),F a (f)]_Q£(t), (10) 

where Q^{t) is a Heisenberg bath operator including the influence of the dynamical subsys- 
tem. Substituting the expression of Eq. (|3|) for this operator, and taking account of the 
fact that the operator of the dynamical subsystem is commutative with Q^{t) at the same 
moment of time (only), we obtain 

A(t) = l[A(t), jy ]_ - l[Q° a (t),Y a (t)} + - l -J dh^itMWait)^^)]^ (11) 

where 

Y a (t) = ±_[A(t),F a (t)U. 

Considering Y a (t) as a function of the unperturbed bath variables and employing the quan- 
tum analog of the Furutsu-Novikov theorem || , we find that 



■\Ql(t),Y a (t)} + ) = J dhM^hX-^^), (12) 

and eliminate the appearance of the bath variables using functional derivatives as described 
by Efremov and Smirnov || 

' <))J1) ) = ( l 7:[Y«(t),Fp(t 1 )]j V (t-t 1 ). (13) 



Introducing the fluctuation source, £(£), 



at) = -\[Q° a (t),Y a (t)} + + J dhM^fatJjp^Fpih)]-, (14) 
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with zero mean value ((£(£)) = 0), we find that the equation for the arbitrary operator A(t) 

of the dynamical subsystem has the form 

1 1 f°° 

A(t) = -[A(t),H ]- --J d*i^(t,ti)[r Q (0,^(*i)] + - 

-J dt 1 M O 0(t,t 1 )jp a (t),Fp(t 1 )]- + Z(t). (15) 

In the case of weak coupling between the dynamical subsystem and the bath, the correlator 
of the fluctuation force is given by 

\ [£(*), £(*!)]+) = M aP {tM) (^ a (i),W] + ) + i«Mi) (^(t),^(ti)]-) , (16) 



where 



#a/?(Mi) = <^[Q°(*),Q?(*i)]-> = ^W(Mi) -</WM))- (17) 



III. FORMULATION OF SPIN DYNAMICS 

We start from the model Hamiltonian 

if = Horbitai + H S pi n + U e - p h + V e-i + UdP + (18) 

describing an electron with spin in the presence of a magnetic field directed along the z-axis, 
where the magnetic field and its vector potential are given by 

B = (0,0,5), A=(-^;^;0). (19) 

The first term in the Hamiltonian fll8"|) is responsible for kinetic electron orbital motion, 

Horbitai = ^ + ^ + 11^, ( 20) 

where the velocity component operators of the electron in a magnetic field can be written 

as 

1 / muj r y\ 1 / muj r x\ p z , . 

m 
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The next term, describing spin motion in the presence of a magnetic field (Zeeman term), 
is given by 

1 . fa^B 

H spin = t;9Vb {cr B) = -^-cr 2 , (22) 

where g is the crystal ^-factor (in particular, g = —0.44 for GaAs), fi B = |e| h/2m c is the 
Bohr magneton (m is the mass of a free electron, the effective mass, m, is 0.067mo for 
GaAs), uj b = gfi B B/h is the frequency of the spin precession induced by the magnetic field, 
and 1? = (a x , a y , <j z ) is the standard set of Pauli matrices. 
The free phonon Hamiltonian is given by 

H B = Y J ^K^ + \), (23) 

k 

where Hlv^ is the phonon energy and and c\ are the creation and annihilation operators of 
a phonon in the k-mode, respectively. The interaction of an electron with random impurities 
is described by the potential 

U^ i (r) = ~^ t U 1t e 1kr , (24) 

k 

where are the spatial Fourier components of the impurity potential, and L 3 denotes the 
volume of the crystal. The electron-phonon interaction is given by 

U e . ph (r,t) = -j^ £<C(6k(f) " ^k(0K kr , (25) 

k 

where ( presents the strength of the electron-phonon coupling. As mentioned above, the 
principal mechanism of spin-orbit interaction is the D'yakonov-Perel' term, Udp- In semicon- 
ductor crystals lacking inversion symmetry, the effective mass Hamiltonian includes terms 
cubic in electron quasimomentum |K| . The presence of these terms induces a random mag- 
netic field O given by 

_ alt , , 

ft = ; r=, 26 

where 



k x = m 3 V x (V; - V 2 ) , x y = m 3 V y {V 2 - V 2 ) , = m 3 V z {V x 2 - V 2 ) , (27) 

Eg is the energy gap, and a is a coefficient representing the strength of the spin-orbit coupling 
(a ~ 0.07 in GaAs). The random magnetic field induced by the orbital motion interacts 
with electron spin through the interaction Hamiltonian 

U DP = \ (o^n) (28) 

It is evident from the Hamiltonian fll8|) that there is no direct coupling of electron spin 
degrees of freedom to phonons and impurities. Therefore, to analyze spin relaxation dy- 
namics, we consider the orbital degrees of freedom as an intermediary that (a) transfers 
fluctuations from the phonon/impurity bath to electron spin, and (b) transfers energy from 
spin degrees of freedom to heat bath. This two-step process is illustrated in Figure 2. Cor- 
respondingly, we address the analysis in two stages, and in the first stage we determine the 
statistical characteristics of orbital motion dictated by the phonon/impurity bath, neglect- 
ing the presence of the electron spin. The pertinent first-stage Hamiltonian, is given 
by 

= HorUtal + Etint + Hph) (29) 

where the first-stage interaction Hamiltonian, H^ t , has the form 

H& = U e . ph + U e ^ = -L~ 3 ' 2 (<C(6k(*) - btjt)) + U k ) e lkr . (30) 

k 

Comparing this expression to Eq. (0), we take the summation index a of Eq. @, as the 
mode number, k, and the function L _3//2 e* kr plays the role of F a (t) of Eq. (0), whereas the 
bath variable, Q a {t), is (i((b k (t) - b± k (t)) + U k ). 

In the second stage of analysis we employ orbital motion as an "effective heat bath" 
having the characteristics determined in the first stage. The second-stage Hamiltonian, 
can be written as 

HM = H spm + + H bath , (31) 



where H\ nt = Udp and Ht, at h may be taken as the remainder of the Hamiltonian (18). To 
compare Udp to Eq. (fj), we rewrite Udp as 



U D p = -cr x Q x (t) - a y Q y (t) - a z Q z (t), (32) 



where the "effective heat bath" variables are given by 



Q*(t) = -^=V x {t) (V y 2 (t) - V*(t)) , (33) 

l\/ ZEg 

3/2 

Q,W = -^=nW (v z 2 (t)-v x 2 (t)), 

Q#) = -%^V,(t) (V x \t)-V*(t)). 

In this, there is summation over coordinate projections and the "dynamic subsystem" vari- 
ables are the electron spin projections. 
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IV. STAGE 1: ORBITAL DYNAMICS. 

In this step we employ TOQS to determine the statistical characteristics of electron 
orbital motion described by the operator equations || 

(it + 7o ) V *® + (Wc + 5) Vyit) = Ut) ' (34) 

J t + 70 ) Vy(t) - (U C - 5) V x (t) = £y(t), 



and 



i + iz)v z (t) = &(*), 



dt 

where V x (t), V y (t), V z (t) are electron velocity operator components, and u> c = \e \ B/mc is the 
cyclotron frequency. The electron-bath interaction determines the relaxation rates, 7o,7z, 
the frequency shift, 5, and the fluctuation sources, £x(t),£ y (t),£ z (t), involved in Eq. (j34|) , 
||. The Fourier transforms of the velocity correlation functions are given by 

\ [KM; v.] + ) = (l \y v (u); v v ] + ) = ^ (- ^— - + . - * 2 ) , (35) 

2 I \2 I 2 \(u)-u c ) +7o (^ + ^c) +7o/ 
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(I [KM; V v ] + ) = -(\ [V v (u); V x } + ) = ^1 (- - - - 1 - ) , 

\ 2 I \ l I l% \{u-u c ) +7o {uj + uj c ) + 75/ 

(36) 

and 

(37) 

where 

(38) 

and 

KM = j dit-^e^^ (^-[Ut)M^)\X (39) 

To be specific we take account of the contributions of polar optical phonons, deforma- 
tional acoustic phonons and random impurities to electron orbital dynamics. In a first 
approximation in their (weak) coupling strengths to orbital motion, these contributions can 
be treated as mutually independent and the overall values of the damping rates and the 
time-Fourier transforms of the correlation functions are given by (superscripts OP and AP 
denote the contributions of optical and acoustic phonons, respectively, and I denotes that 
of random impurities) : 

7o = 7o° P + 7o 7 + 7o AP , (40) 

7,=7r+7 2 7 +7r, 

K ± (u) = K° p {u) + K{(u) + Kf(co), 
K z {u) = K° p {u) + Ki{u) + Kf{u\ 

The microscopic expressions of all these quantities are determined using the method of Ref. 
|J, and are given below: 

The response and correlation functions for polar optical phonons are 
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op, \ 47rfi e 2 OP hAitVL Q e 2 

V^k ( T ) = sm(Q r)?7(T), (r) 



M7 



cos (finr) coth(— — — ), 



(41) 



resulting in 



7o 



op _ 1 



2 /-+oo 

/ dk ± k 3 ± 
Jo Jo 



2v / 27r me* J 



r c k z ) op 



k 2 



& UF (k±,k z ) 



and 



7. 



op 1 ^o e 



2 /-+oo /• 

/ dk±k± 
Jo Jo 



2tt me* J 



r c (fcj., k z ) , 2^op 



k 2 



k 2 z G OP (k ± ,k z ), 



where 



G OP (k ± , k z ) = \ (coth(^L) + 1 W + fi ) exp 



-^(cj fc + fio) 2 r c 2 (/c_L,A; 2 ) 



+ (42) 



+ ( coth (^~^) - 1 ) (^fc - fi o) exp 



~ (u k - tt ) 2 r 2 (k ± ,k z ) 



Furthermore, we obtain 
and 



4V/27T m 2 e* jo 



ft 



(43) 



2\/2¥ m 2 e* 7 



p+oo P+OO 

/ dk±k± / 
Jo Jo 



T c (k±, k z ) j OirrOP 



k 2 



k 2 z K UF (k ± ,k z ,oo), (44) 



where 



K op (A; ± ,A; 2 ,cj) 



coth (^^) + 1 ) ( ex P 



--(uj + uj k + tt ) r c (k ± ,k z ) 



+ (45) 



+ exp 



1 



--(u-u k -n ) T c (k±,k z ) 



+ 



+ ( coth (|^) " 1 ) ( exp 



--{u + u k -Vto) r c (k ± ,k z ) 



+ 



+ exp 



~ (uj - u k + tt ) 2 r 2 (k ± , k z ) 



Qo is the optical phonon frequency, 1/e* = l/e^ — 1/eo, (foe and eo are the hf and static 
permittivities of the crystal, respectively). 

The response and correlation functions of the acoustic phonons are given by 
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AP( T ) = _ZA gin ^ ukr) ^ M AP (r) = \— cos (ukr) coth( 7 ^), 
pu 2 pu 2.Kb± 

which provide the contribution to the damping rate as 

1 D 2 r+co 



7o AP 



87TV27T mpu 



(46) 



r+co r+co 

/ dk L k\ / dk z kT^(k ± ,k z )G AP (k±,k z ) (47) 
Jo Jo 



and 



7, 



AP 



1 D 



2 r+oo r+co 

/ rfA; ± A; ± / dk z kk 2 z T*(k ± , k z )G AP (k ± , k z ) , 
Jo Jo 



(48) 



where 



fat A; 



G k z ) = { ( cotM^ry ) + 1 ) K + exp 

fa A; 



Accordingly, 
and 

where 

K ap (A; ± ,A; 2 ,cj) 



(fat A; \ 
coth(^-^) - 1 J (w fc - ufc) exp 



l ///;- /■ x „ 



-- (w fc + uk) 2 r 2 (k ± , k z ) 
~ {uj k -ukf T 2 (k_L,k z ) 



+ (49) 



167rV27rm z pw J 



r+co 

k\ I dk z kr c (k±, k z )K AP (k±, k z , u) , 
Jo 



I fcr)2 /-+00 /-+oo 

/ dfe^ / dk z kk 2 r c (k ± , k z )K AP (k ± , k z ,u) , (51) 
Jo Jo 



87rV27r m^pw J 



c °th(^^) + 1 ) ( exp 



--(c^ + u; fe + wA;) T c (k ± ,k z ) 



+ exp 



--(w-Wfc-ufc) T c (k ± ,k z ) 



+ 



+ ( coth (^) " 1 ) ( ex P 



(a; + t^fc - wA;) 2 t c 2 (/c_l, A^) 



+ 



+ exp 



-- (u - u) k + ukf r 2 (k ± , k g ) 



(50) 



+ (52) 



Here, D is the deformation potential , p is the crystal density, and u is the sound velocity. 

For the case of static random impurities, there is not any response function, due to the 
lack of dynamics. The correlation function is given by 
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Ml 



2e 4 n* 



(53) 



ire* (P + r - 2 ) 2 ' 

where r is the screening radius, n* = J2 a n a^ay n a is the impurity concentration for species 
a, and Z a is their charge number. The relaxation rates associated with impurities have the 
forms 

i _ 1 e 4 n* 

and 



dk i fci 



00 ^ k 2 r?(k ± ,k z ) 
dk z q- exp 

(F + r - 2 ) 2 



uJ 2 k T 2 (k ±} k z ) 



(54) 



7f 



4„* Z'+oo 



en 



7r 2 \/27rm 2 eo Jo 



dk\k\ 



+00 



kzk 2r c{k± } k z ) 



2\ 2 



(fc 2 + r " 2 ) 



exp 



u 2 k r 2 (k ± ,k z ) 



(55) 



Finally, the contributions to the Fourier transforms of the fluctuation force correlators can 
be written as 



exp 



dk] k. 



2\2 



{k 2 + r - 2 ) 



(56) 



-~(u + u k ) T c (k ± ,k z ) 



+ exp 



■-(o;-a; fc ) T c (k ±} k z ) 



and 



e 4 n? 



7r 2 v / 27r m2 eo 



exp 



dk\k\ 



-~ (^ + ^ fe ) 2 r 2 (A;_ L ,/i; z ) 



2\ 2 



(P + r - 2 ) 



(57) 



exp 



-^(w-w A; ) 2 r 2 (A; ± ,^) 



In all these formulae we have used the notation k = \Jk\-\- k 2 , u k = hk 2 /2m, 
T~c 2 (k_i_,k z ) = k\ (V 2 ) + k 2 z (V 2 ). (V 2 ) and (V 2 ) can be determined self-consistently using 
(V 2 ) = K ± (uj c )/2-f and (V 2 ) = K z (0)/2^ z @. All these results are obtained for relatively 
weak magnetic fields, u c r c << 1. However, the energy shift, S, due to electron-bath coupling 
is even smaller, 5 « u> c , and will be neglected. 



V. STAGE 2: SPIN DYNAMICS. 



In this section we examine the spin relaxation process in the presence of the "effective 
bath" with variables given by Eq. fl3"Bp. Their response functions and correlation functions 
have the form (i, j = x, y, z) 
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AMMO = (^QiW'Qi^U) 



(58) 



(59) 



It is evident that these functions are of the sixth order of the electron velocity projections. 
Using the quantum analog of the Furutsu-Novikov theorem for Gaussian variables they 
can be represented in terms of pair correlators. The resulting expressions for the spectral 
functions are given by 

Sij (u)= fd(t- h) e^-^Mij (t, h) , (60) 



S„(u) = S yy (u) = A 2 (- [V x (u),V x ] + j [{(V z 2 ) - (V x 2 )Y + A(V x V y ) (V x V y ) 
2 f duj x f du 2 _ 

+ A / — / —^{UJi,UJ2,CO-COl-UJ2)- 

2 /i [VxM, V x } + \ /i [V x {u 2 ), V x } + \ /- [V x (uj -Ui- w 2 ), 14]. 
2 /i [y^), ^]+) (5 [^(u; 2 ), T4] + ^ ^ [T,(w - wi - u*), VJ_ 
-4 /l [KM, /i [T4(u; 2 ), V v ] + \ /i [KV - W! - a*), V y ]_ 



+ 



(61) 



S xy {uj) 



and 



-^M = a 2 / - \y y (v), v x ] + ) [«k 2 > - (V x 2 )) 2 + 4 (V x V y ) (V x V y ) 

2 /" dwx f cL> 2 



+ 



2 L-« 



[^(^),K.]_ 



P4(u; - wi - w 2 ), ) + 



+2(i[r y (a^) > ^] + \/i[K(w a ) > vj. 



- [k(^ -w x -w 2 ),y*]. 



SUM = A 



, duji f duj 2 „ 

s (a;i,a; 2 ,a; — cji — a; 2 J 



2tt / 2tt 



■4(-[K,(w 1 ),yj. 



[K(w 2 ),Kl 



[14 (a; - wi - co 2 ), V x ], ) - 



(62) 



(63) 



- [V x (uJ 2 ),Vy] + \ /- [14 (W - Wl - W2), 14]. 
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where we defined an auxiliary function 

5 (ui, LO2, W3) = 1 + tanh ( — =7 J tanh (—^ ) + 



2T \2T 



+ tanh tanh (^) + tanh (^) tanh , (64) 

which appears due to the replacement of the velocity commutators by their anticommu- 
tators in accordance with Eq. ([|). In Eqs. (|6l| - |63D we use the shortened notation 

(V?> = (U v Mv j (t)) + ) = St{WM,Vi\ + ), ^d (v x v y ) = (im),v y (t)} + ) = 

I irr (2 [^( t<; )' + ) -Employing a second application of the TOQS method, we obtain equa- 
tions for the average spin projections as follows: 

j t (a x (t)) = -T xx (a x (t)) - (u B + V xy ) (a y (t)) + T xz (a x (t)) + T° x , (65) 
j t M*)> = (wb + T,,) (a x (t)) - T yy (a y (t)) + Y yz (a z (t)) + T° y , 
j t (a z (t)) = Mt)) + T zy (a y (t)) - T zz (a z (t)) + T° z , 

where the coefficients are given by 

T xx = ^ J dti {M yy (t, h) cos [u B [t - t x )} + M yx (t , t x ) sin [u B (t - h)] + M zz (t, h)} , (66) 
4 /"* 

T xy = J dti {M yy (t, ti) sin [wb (t - ti)] - (t, tj) cos [u B (t-ti)]}, 
4 /•* 

r xz = ^ y (*' *0 cos (* ~~ _ M *v fa ti) sin [u B {t - h)]} , 

2 f* 

r° = - J dti Wzx (t, h) sin [uj b {t - h)) + (p zy (t, tt) cos [u B {t - h)) - (p yz (t, tt)} , 



r = 

1 yy 

r =± 

2/2 



T yx = ^ j dh {M xx (t, h) sin [u B {t - t t )\ + M xy (t, h) cos [u B (t - h)]} , (67) 
4 /•* 

— / dt x {M xx {t, t x ) cos [u B (t - h)} - M xy {t, h) sin (t - t x )] + 
4 r* 



r»- 2 

and 



o —00 

^ y d*! {M^ (*, ti) sin [u B (t - h)] + M zy (t, h) cos [u B (t - tj}} , 

2 r* 

- J dt x {ip zy (t, ti) sin [u B (t - ti)} - ip zx (t, ti) cos [w fl (t - ti)] + ^ (t, ti)} , 
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T zx = A I dhM.-At-ii). (68) 



OO 



4 

r 2y = ^ y dtiM va (t,ti) , 

4 

r 22 = ^ y d*i {(M ra (t, h) + M ra (f , *x)) cos [uj B (t + 

+ {My X (t, h) - M xy {t, h)) sin [u B (t-h)]}, 
r° = ^ y dti {(w (*, *i) + ^ (*, *i)) cos [lo b (t - A)} - 
- (<Pxx (t, h) + ipy y (f , tx)) sin [u B (t - ti)]} • 

Only six of the above twelve coefficients, T^, do not vanish, and we may identify relaxation 
times in terms of the nonvanishing T^- as follows: 

^xx — Fyy — T^7i ^ zz = 7^7) (69) 
J-2 1 1 

where T 1? T 2 are the longitudinal and transverse relaxation times, responsible for the relax- 
ation of the magnetic moment and decoherence respectively. Furthermore, direct calculation 
shows that 

T xy = T yx = 5 S , T° z = a° z = - tanh (hw B /2T) , (70) 

where 5 S may be identified as a frequency shift due to "effective bath" fluctuations, and a z 
is the equilibrium population difference between spin up and spin down states. With these 
calculated results and identifications, the spin equations take the usual Bloch form, 

j t <**(*)> = " ("b + &.) Mt)) , (71) 

±(a y (t))=(u> B + 5 a )(a x (t))-£^, 

d Irr ^ ~ M*)) 

di Mt)) = T\ • 

In terms of the spectral functions, Sij(cu), the relaxation rates (inverse relaxation times), 
may be expressed as follows: 

1 4 

7fT = £2 ( S «N + iS xyM) , (72) 
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and 



T 2 h 2 



(S xx (ujb) + iS X y(uj B ) + S zz (0)) 



(73) 



Although Eqs. (|72|), (^) are relatively simple and yield some interesting qualitative conclu- 
sions, a quantitative analysis of them is difficult, because of the complexity of the expressions 
for the spectral functions, Sij(uj). This complexity is relieved under the prevailing assump- 
tion of weak coupling, which permits the replacement of the Lorenzians involved in the 
integrands of Eqs. (|B~T1 - |6"3"|) by 

1 e 



A(x;e) 



5(x), when e — » 0, 



7T x 2 + e 2 

where S(x) is the Dirac delta-function, and, consequently, 

dxf(x)A(x; ea)A(x - y; eb) = ^A(y; e (a + b)) (f(y) + /(0)) + 

+lA(y;e(o-6))(/(y)-/(0)) 



(74) 



(75) 



(In this, we note that in Eqs. fl35]),(|3~6D and (|37|), we have K±(u), K z (u), jo, j z all propor- 
tional to coupling strength ~ e). Having employed Eqs. ([7|), (|75|) in the integrals of Eqs. 
representing S xx (u), S xy (u), and S zz (u>), we obtain the results of integration as: 



^(^ VrUlV,..) ((V z 2 )-(V x 2 )) 2 + A(V x V y )(V x V y ) 



SE g \2 

2 S 

am 



+ 



(76) 



8E a 



[i2i(w;w c ) + Y x {u;u c )] 



S X y{w) 



a 2 m 3 1 1 



9 3 

am 



8£„ 



- [V X (UJ), Vy} + ) {(V z 2 ) - (V X 2 )Y + 4 (V X Vy) (V X Vy) 



[R 2 (uj;uj c ) + Y 2 {uj\u c )} 



(77) 



2 3 

am 7r 



-K_l(o; c 



{A(cj - 2w c ; 270 + 7 Z ) (A" z (0)A" ± (a; - u c )E(u - u c ; u c ; 0) + 
+K z (u - 2a; c )if_i_(u; c )H(u; - 2w c ; u c ; lo c )) + 



(78) 
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+A(u - 2u c ; 270 - 7z ) (K z (0)K±(u - lo c )E(uj - u c ; u c ; 0)- 

-K z (u - 2uj c )K 1 _(uj c )Q(uj - 2u c ; uj c ; uj c )) + 

+A(cu + 2oo c ] 2 7o + 7*) (7^(0)7^ (cj + cj c )S(cj + cj c ; -cu c ; 0) + 

+i^ z (cu + 2u c )K ± (u c )E(u + 2cj c ; -cj c ; -u c )) + 

+A(cj + 2cj c ; 2 7o - 7 Z ) (^(0)ir ± (c; + cj c )S(cj + cj c ; -cu c ; 0)- 

-K z {uj + 2<^ c )lf_ L (u; c )S(u; + 2cj c ; -u c ; -u c ))} , 

where we have introduced the notation 

TT 2 

i2i(w,w c ) = —^7 (K_l(cj c )) 
I670 

{3A(w - 3cj c ; 3 7o )-ft'±(cj - 2cj c )S(cu - 2w c ; uj c ) + 
+A(u - u c ; 3jo)K±(u - 2u c )E(u - 2u c ; io c \ io c ) + 

+2A(cj - uj c ; 70) (K ± (u)E(u; co c ; -uj c ) - K ± (uj - 2u c )E(u - 2u c ; uj c ; uj c )) + 
+2A(cj + u c ; 70) (K± (u) E(u; lo c ; -uj c ) - K ± (u + 2uj c )E(uj + 2u c ; -u c ; -u c )) + 
+A{oj + oj c ] 37o)i^±(o; + 2l} c )E(u + 2oj c ; —u c ; —uj c ) + 
+3A{uj + 3w c ; 3'y )K±(u! + 2u c )E(u + 2uj c \ —ou c ; — uj c )} , 

R 2 {u,uj c ) = -i— m (K±(cj c )) 2 
16 7 o 

{3A(oj — 3u c ; 3 r y )K±(uj — 2uj c )E{uj — 2uj c ; u c ; u c ) — 
-A(uj - u c ; 3^ )K ± (uj - 2uj c )E(uj - 2uj c ; uj c ; uj c ) - 

-2A(cj - u c ; 70) (K ± (u)E(u; uj c ; -uj c ) - K±(u - 2uj c )E(uj - 2uj c ; uj c ; uj c )) + 
+2A(u + u c ; 70) (K ± (u)E(uj] u c ; -u c ) - K±(u + 2u c )E(u + 2u c ; -u c ; -o> c )) + 
+A(u + uj c \ 37o)i^j_(o; + 2cj c )S(cj + 2u c ; —oj c ; —u c ) — 
-3A{uj + 3uj c ; 3-f Q )K 1 _(uj + 2uj c )E(uj + 2u c ; -uj c ; -uj c )} , 

{A(u - u c - 2 lz + 70) (K ± (u c )K z (u - u c )E{u - u c ; u c ; 0) + K ± (u)K z (0)) + 
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+A(u - u c ; 2^ z - 70) (K±(uj c )K z (uj - uj c )E(u - u c ; u c ; 0) - Kj_(u)K z (0)) + 
+A(u + u c ; 27 2 + 70) (k±(-u) c )K z (u + uj c )E(uj + u c ; -u c ; 0) + K xx (u)K zz (0) J + 
+A(w + u c ; 2 lz - 7o ) (ir ± (-w c )^(cj + lu c )E(uj + u c ; -u c ; 0) - K ± (u) K 2 (0))} , 



and 



Z7T 



^7 2 7o 



i^ 2 (0) 



(79) 



{A(cj - u c - 2 7z + 7o ) (if ± (u; c )i^(u; - uj c )E(u - u c ; u c ; 0) + K ± (u)K z (0)) + 
+A(u - u c ; 27 2 - 7 ) (K ± (u e )K x (u - uj c )E(u - w c ; u c ; 0) - J fT_ L (u;)^(0)) - 
-A(oj + lu c ; 2 lz + 7o ) (K ± (-lu c )K z (lu + cu c )H(cu + uj c ; -u e ; 0) + iT ± (cj)^(0)) - 
-A(u + cu c ; 2 7z - 7o ) (i^-u^i^u; + cj c )H(w + cu c ; -w c ; 0) - K ± (u)K z (0))} . 

These results are very general, providing explicit expressions for the spin relaxation times in 
terms of the material constants and coupling strengths. The magnetic field has been treated 
as relatively weak {uj c t c « 1), but our estimates show that this range of fields includes 
those used in most experiments. 

To make contact with earlier theories that did not take account of the magnetic field 
and used a phenomenological momentum scattering time |K|, we put u c — and ujb = in 



Eqs. (|76|-[79|), which immediately yields 



where 



C l \ O ( \ Q l \ « 2 ™ 3 7T [if (0)f Jjrl v., v 

S xx {u) = S xx {uj) = S zz {u) = — 3 K [uj)A(uj; 3a~f ) 

°&g 7o 



K ± (u) = K z (uj) = K(u), 70 = 72 



(80) 



(81) 



and 



S X y(uj) = 0. 

In this zero field limit the relaxation times are equal and determined by 

Ti T 2 h 2 



(82) 



(83) 
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In equilibrium, the rate of velocity fluctuations is defined only by temperature, and we 
obtain the result 

£(o) kjr 

V •' 2 7 „ m ' 1 ' 

where T is the Kelvin temperature. This yields the common relaxation rate as 

1 4 a 2 1 a 
Ti >2 3 70 

which may be brought into coincidence with the expression usually used for the DP mecha- 
nism f7|,j8|, |10l] by the change of notation qr p — > g — , with r p a phenomenological momentum 
scattering time. It should be noted that the choice of coefficient q is unclear, with various 
authors assigning values in the range q = 0.8 — 2.7, with q ph 2.7 for scattering by defor- 
mational acoustic phonons and q ~ 0.8 for scattering by polar optical phonons [§H[U|. On 
the other hand, our microscopic analysis provides a clearly defined value with account of all 
scattering mechanisms. Moreover, we provide results for finite magnetic field derived on a 



microscopic basis, rather than relying on the ansatz of Ref. ||10|1 , which has T\(B) = T%(B), 

1 _ 1 1 

T^B) ~ T li2 (fl = 0)l + (u B T p f 1 } 

In contradistinction to this, our analysis in the presence of a magnetic field shows that the 
transverse and longitudinal spin relaxation times are not equal (Eqs. (|72|,[73D). 



VI. SUMMARY 

In Summary, we have analyzed electron spin relaxation dynamics and decoherence in bulk 
semiconductors with an applied external magnetic field. In the absence of direct coupling of 
spin degrees of freedom to the phonon/impurity heat bath, we have used the orbital motion 
as an intermediary, transferring fluctuations from the bath to electron spin and, transmitting 
the dissipated energy from spin to the bath. To accomplish this, the two-stage procedure of 
solution has been employed: In the first stage, the fluctuation characteristics of the orbital 
motion were determined and the orbital degrees of freedom were used in the second stage as 

20 



an "effective heat bath" for spin dynamics. In both stages of the analysis, we employed the 
method of the theory of open quantum systems, obtaining a set of Bloch equations having 
longitudinal and transverse spin relaxation times, T\ and T 2 , determined on a microscopic 
basis. Our results provide analytic definition of the phenomenological parameters employed 
in earlier zero-field theories, in terms of material constants and coupling strengths. Moreover, 
this analysis of electron spin dynamics in finite magnetic field yields explicit formulae for 
the longitudinal and transverse relaxation times. 

Aknowledgement 

V.I. P., L.G.M. and N.J.M.H. gratefully acknowledge support from the Department of 
Defense, DAAD 19-01-1-0592. 



21 



FIGURES 

FIG. 1. Schematic of "dynamical subsystem" - "heat bath" interaction. 
FIG. 2. Schematic of the two-stage procedure for spin relaxation analyses 
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